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We present a measurement of the time-dependent CP- violating asymmetry in _B° — >■ J/tp (f) decays, 
using data collected with the LHCb detector at the LHC. The decay time distribution of — >■ J/^j) 4> 
is characterized by the decay widths Fh and Fl of the heavy and light mass eigenstates of the B^-B^^ 
system and by a CP- violating phase (j)s. In a sample of about 8500 5° — >■ J/ip events isolated from 
0.37 fb~^ of pp collisions at = 7TeV we measure <f)s = 0.15 ± 0.18 (stat) ± 0.06 (syst) rad. We 
also find an average B° decay width Vs = (Fl + Fh)/2 = 0.657 ± 0.009 (stat) ± 0.008 (syst) ps"^ 
and a decay width difference AF^ = Fl - Fh = 0.123 ± 0.029 (stat) ± 0.011 (syst) ps~\ Our 
measurement is insensitive to the transformation {<j)s, AFs) (tt — (ps, — AFs). 

To be submitted to Physical Review Letters 



In the Standard Model (SM) CP violation arises through 
a single phase in the CKM quark mixing matrix [1 . In 
neutral B meson decays to a final state which is ac- 
cessible to both B and B mesons, the interference be- 
tween the amplitude for the direct decay and the ampli- 
tude for decay after oscillation, leads to a time-dependent 
CP-violating asymmetry between the decay time distri- 
butions of B and B mesons. The decay B^ J/4'4> 
allows the measurement of such an asymmetry, which 
can be expressed in terms of the decay width differ- 
ence of the heavy (H) and light (L) _B" mass eigen- 
states AFs = Fl — Fh and a single phase (t)s 2 . In 
the SM, the decay width difference is AF™ = 0.087 ± 
0.021 ps~^ [3 , while the phase is predicted to be small, 
(t^f^ = -2arg(-14,y;,/KsV;;) = -0.036±0.002 rad 0]. 
This value ignores a possible contribution from sub- 
leading decay amplitudes [5] . Contributions from physics 
beyond the SM could lead to much larger values of i^g [B] . 

In this Letter we present measurements of (ps , AF^ and 
the average decay width Tg = (Fl + Fh)/2. Previous 
measurements of these quantities have been reported by 
the CDF and D0 collaborations [7^. We use an integrated 
luminosity of 0.37 fb^^ of pp colhsion data recorded at a 
centre-of-mass energy -^s — 7TeV by the LHCb experi- 
ment during the first half of 2011. The LHCb detector 
is a forward spectrometer at the Large Hadron Collider 
and is described in detail in Ref. [8]. 

We look for 5° — ^ J/i^4> candidates in decays to 
J/ij} — > ijT and cf) — > K^K^ . Events are selected by 
a trigger system consisting of a hardware trigger, which 
selects muon or hadron candidates with high transverse 
momentum with respect to the beam direction (px), fol- 
lowed by a two stage software trigger. In the first stage 
a simplified event reconstruction is applied. Events are 
required to either have two well-identified muons with 
invariant mass above 2.7 GeV, or at least one muon or 
one high-pT track with a large impact parameter to any 
primary vertex. In the second stage a full event recon- 
struction is performed and only events with a muon can- 
didate pair with invariant mass within 120 MeV of the 
nominal J/ip mass [S] are retained. We adopt units such 
that c = 1 and h = 1. 

For the final event selection muon candidates are re- 
quired to have p^ > 0.5 GeV. J/i/' candidates are cre- 



ated from pairs of oppositely charged muons that have 
a common vertex and an invariant mass in the range 
3030 — 3150 MeV. The latter corresponds to about eight 
times the /i"*"^" invariant mass resolution and covers part 
of the J/'ip radiative tail. The selection requires two 
oppositely charged particles that are identified as kaons, 
form a common vertex and have an invariant mass within 
±12 MeV of the nominal mass [9 . The of the 
candidate is required to exceed 1 GeV. The mass window 
covers approximately 90% of the (p — )■ K~^K~ lineshape. 

We select 5° candidates from combinations of a J/ip 
and a (j) with invariant mass m the range 5200 — 
5550 MeV. The latter is computed with the invariant 
mass of the /i^/i^ pair constrained to the nominal J ftp 
mass. The decay time t of the B'^ is obtained from a ver- 
tex fit that constrains the 5° — )■ jX^ fjT K~ candidate 
to originate from the primary vertex [10] . The of the 
fit, which has 7 degrees of freedom, is required to be less 
than 35. In the small fraction of events with more than 
one candidate, only the candidate with the smallest 
is kept. B^ candidates are required to have a decay time 
within the range 0.3 < t < 14.0 ps. Applying a lower 
bound on the decay time suppresses a large fraction of 
the prompt combinatorial background whilst having a 
small effect on the sensitivity to (ps- From a fit to the 
nriB distribution, shown in Fig. [T] we extract a signal of 
8492 ± 97 events. 

The B° — >■ J/ip(j) fj,^ decay proceeds via 

two intermediate spin-1 particles (i.e. with the K^K^ 
pair in a P-wave). The final state can be CP-even or 
CP-odd depending upon the relative orbital angular mo- 
mentum between the J/tp and the (p. The same final 
state can also be produced with K~ pairs with zero 
relative orbital angular momentum (S-wave) [11] . This 
S-wave final state is CP-odd. In order to measure (ps 
it is necessary to disentangle the CP-even and CP-odd 
components. This is achieved by analysing the distribu- 
tion of the reconstructed decay angles — {9,ip,ip) in 
the transversity basis [HI |T3] . In the J/ip rest frame we 
define a right-handed coordinate system such that the x 
axis is parallel to the direction of the (p momentum and 
the z axis is parallel to the cross-product of the K~ and 
momenta. In this frame 9 and ip are the azimuthal 
and polar angles of the The angle ip is the angle 
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FIG. 1. Invariant mass distribution for — >■ fi^ fj,^ K'^ K~ 
candidates with the mass of the n'^ ^~ pair constrained to 
the nominal J/tji mass. Curves for fitted contributions from 
signal (dashed), background (dotted) and their sum (solid) 
are overlaid. 



between the K~ momentum and the J/i}) momentum in 
the rest frame of the 

We perform an unbinned maximum likelihood fit to the 
invariant mass m^, the decay time t, and the three decay 
angles 51. The probability density function (PDF) used 
in the fit consists of signal and background components 
which include detector resolution and acceptance effects. 
The PDFs are factorised into separate components for 
the mass and for the remaining observables. 

The signal distribution is described by two Gaus- 
sian functions with a common mean. The mean and 
width of the narrow Gaussian are fit parameters. The 



fraction of the second Gaussian and its width relative to 
the narrow Gaussian are fixed to values obtained from 
simulated events. The tub distribution for the combina- 
torial background is described by an exponential func- 
tion with a slope determined by the fit. Possible peaking 
background from decays with similar final states such as 
— > J/tpK*'^ is found to be negligible from studies 
using simulated events. 

The distribution of the signal decay time and angles 
is described by a sum of ten terms, corresponding to the 
four polarization amplitudes and their interference terms. 
Each of these is the product of a time-dependent function 
and an angular function |12) 
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k = l 



The time-dependent functions hk (t) can be written as 

hk{t) = Nke-^^' 
+ak 



[ck cos(A'mst) +dkSm{Amst) 
cosh (^AF^t) -t- bk sinh (jAFji)] . (2) 

where Anis is the oscillation frequency. The coeffi- 
cients Nk and au, ■ ■ ■ ,dk can be expressed in terms of (ps 
and four complex transversity amplitudes Ai at t = 0. 
The label i takes the values {_L,||,0} for the three P- 
wave amplitudes and S for the S-wave amplitude. In the 
fit we parameterize each ^4^(0) by its magnitude squared 
|Ai(0)p and its phase 6i, and adopt the convention 6o — 
and = 1- For a particle produced in a B° 

flavour eigenstate the coefficients in Eq. [2] and the angu- 
lar functions fk {^) are then, see [T31 [H] , given by 



k 




Nk 




bk 


Ck 


dk 


1 


2 cos^ V' (l ~ sin^ cos^ 4>) 


l^o(0)P 


1 




— cos (ps 





sin 0s 


2 


sin^ (l ~ sin^ 9 sin^ 0) 


l^||(0)P 


1 




— cos (ps 





sin (ps 


3 


sin^ V sin^ 6 




1 




COS (ps 





- sin 0s 


4 


— sin^ ip sin 29 sin (f> 


(0)^^(0)1 







— cos((5j_ — ) sin ips 


sin(Jj_ — <5||) 


— cos((5j_ — ) cos 0s 


5 


i^Ssin 2->p sin2 6'sin20 


l^o(OM||(0)| 


cos((5|| — 


So) 


— cos((5|| — So) cos 0s 





cos((5|| — 6o) sin 0s 


6 


i \/2 sin 2tp sin 29 cos (j) 









— cos((5j_ ~ <5o) sin (ps 


sin((5_L - So) 


— cos((5j^ — 6o) cos 0s 


7 


1(1 -sin^ffcos^ (j)) 
^ -\/6 sin ip sin^ 9 sin 2(j) 




1 




cos (ps 





— sin 0s 


8 


l^s(0)A||(0)| 







— sin((5|| — (5s) sin pg 


cos(^ll — 6s) 


— sin((5|| — 6s) cos 0s 


9 


^ -s/G sin ip sin 29 cos (j) 


\Asm±m 


sin((5j^ — 


Ss) 


sm{6± — 6s) cos (ps 





— sin(5j^ — 6s ) sin 0s 
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|\/3cos'0(l — sin^ ^cos^ 0) 


|^s(0)Ao(0)| 







- sin((5o - Ss) sin 0s 


cos((5o - 6s) 


-sin((5o -(5s) cos 0s 



We neglect CP violation in mixing and in the decay 
amplitudes. The differential decay rates for a B^ meson 
produced at time t — are obtained by changing the 
sign of 0s, ^_l(0) and As(0), or, equivalently, the sign 



of Ck and dk in the expressions above. The PDF is in- 
variant under the transformation (0s, AFs, (5|| , (5j_, ^s) ^ 
(tt — 0s, — AFs, — (5|| , TT — (5j^,— 5s) which gives rise to a 
two-fold ambiguity in the results. 
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We have verified that correlations between decay time 
and decay angles in the background are small enough to 
be ignored. Using the data in the rriB sidebands, which 
we define as selected events with niB outside the range 
5311 — 5411 MeV, we determine that the background 
decay time distribution can be modelled by a sum of two 
exponential functions. The lifetime parameters and the 
relative fraction are determined by the fit. The decay 
angle distribution is modelled using a histogram obtained 
from the data in the niB sidebands. The normalisation of 
the background with respect to the signal is determined 
by the fit. 

The measurement of requires knowledge of the 
fiavour of the meson at production. We exploit the 
following flavour specific features of the accompanying 
(non-signal) 6-hadron decay to tag the fiavour: the 
charge of a muon or an electron with large transverse mo- 
mentum produced by semileptonic decays, the charge of 
a kaon from a subsequent charmed hadron decay and the 
momentum-weighted charge of all tracks included in the 
inclusively reconstructed decay vertex. These signatures 
are combined using a neural network to estimate a per- 
event mistag probability, w, which is calibrated with data 
from control channels [15] . The fraction of tagged events 
in the signal sample is £tag = (24.9 ±0.5)%. The dilution 
of the CP asymmetry due to the mistag probability is 
D = 1 — 2lo. The effective dilution in our signal sample 
is D ^ 0.277 ±0.006 (stat)± 0.016 (syst), resuhing in an 
effective tagging efficiency of etag-D^ = (1.91 ± 0.23)%. 
The uncertainty in uj is taken into account by allowing 
calibration parameters described in Ref. |15j to vary in 
the fit with Gaussian constraints given by their estimated 
uncertainties. Both tagged and untagged events are used 
in the fit. The untagged events dominate the sensitivity 
to the lifetimes and amplitudes. 

To account for the decay time resolution, the PDF is 
convolved with a sum of three Gaussian functions with a 
common mean and different widths. Studies on simulated 
data have shown that selected prompt J/ipK'^K~' combi- 
nations have nearly identical resolution to signal events. 
Consequently, we determine the parameters of the res- 
olution model from a fit to the decay time distribution 
of such prompt combinations in the data, after subtract- 
ing non-J/i/; events with the sPlot method ,16j using the 
invariant mass as discriminating variable. The re- 
sulting dilution is equivalent to that of a single Gaussian 
with a width of 50 fs. The uncertainty on the decay time 
resolution is estimated to be 4% by varying the selection 
of events and by comparing in the simulation the reso- 
lutions obtained for prompt combinations and signal 
events. This uncertainty is accounted for by scaling the 
widths of the three Gaussians by a common factor of 
1.00 ± 0.04, which is varied in the fit subject to a Gaus- 
sian constraint. In similar fashion the uncertainty on 
the mixing frequency is taken into account by varying it 
within the constraint imposed by the LHCb measurement 



Ams = 17.63 ± 0.11 (stat) ± 0.02 (syst) ps"! 

The decay time distribution is affected by two accep- 
tance effects. First, the efficiency decreases approxi- 
mately linearly with decay time due to inefficiencies in 
the reconstruction of tracks far from the central axis 
of the detector. This effect is parameterized as e(t) oc 
(1 — j3t) where the factor /? = 0.016 ps~^ is determined 
from simulated events. Second, a fraction of approxi- 
mately 14% of the events has been selected exclusively 
by a trigger path that exploits large impact parameters 
of the decay products, leading to a drop in efficiency at 
small decay times. This effect is described by the empiri- 
cal acceptance function e{t) oc {atY / [1 ± [atY], apphed 
only to these events. The parameters a and c are deter- 
mined in the fit. As a result, the events selected with 
impact parameter cuts do effectively not contribute to 
the measurement of Tg ■ 

The uncertainty on the reconstructed decay angles is 
small and is neglected in the fit. The decay angle accep- 
tance is determined using simulated events. The devia- 
tion from a flat acceptance is due to the LHCb forward 
geometry and selection requirements on the momenta of 
final state particles. The acceptance varies by less than 
5% over the full range for all three angles. 

The results of the fit for the main observables are 
shown in Table ITl The likelihood profile for 5\\ is not 
parabolic and we therefore quote the 68% confidence level 
(CL) range 3.0 < 5\\ < 3.5. The correlation coefficients 
for the statistical uncertainties are p(rs,Ars) = —0.30, 
piTs,(l)s) = 0.12 and p{ATs,(t>s) = -0.08. Figure [5] shows 
the data distribution for decay time and angles with the 
projections of the best fit PDF overlaid. To assess the 
overall agreement of the PDF with the data we calculate 
the goodness of fit based on the point-to-point dissimilar- 
ity test [ISJ. The p- value obtained is 0.68. Figure [s] shows 
the 68%, 90% and 95% CL contours in the AF^-^^ plane. 
These contours are obtained from the likelihood profile 
after including systematic uncertainties, and correspond 
to decreases in the natural logarithm of the likelihood, 
with respect to its maximum, of 1.15, 2.30 and 3.00 re- 
spectively. 

TABLE I. Fit results for the solution with AF^ > with 
statistical and systematic uncertainties. 



parameter 


value 


Cstat. 




Ls [ps~^] 


0.657 


0.009 


0.008 


Ar, [ps-1] 


0.123 


0.029 


0.011 




0.237 


0.015 


0.012 


l^o(0)P 


0.497 


0.013 


0.030 


\Asm' 


0.042 


0.015 


0.018 


5± [rad] 


2.95 


0.37 


0.12 


5s [rad] 


2.98 


0.36 


0.12 


(ps [rad] 


0.15 


0.18 


0.06 



3 




' COS V (p [rad] 



FIG. 2. Projections for the decay time and transversity angle 
distributions for events with ms in a ± 20 MeV range around 
the B° mass. The points are the data. The dashed, dotted 
and solid lines represent the fitted contributions from signal, 
background and their sum. The remaining curves correspond 
to different contributions to the signal, namely the CP-even 
P-wave (dashed with single dot), the CP-odd P-wave (dashed 
with double dot) and the S-wave (dashed with triple dot). 

The sensitivity to (j)s stems mainly from its appear- 
ance as the amplitude of the sin(Amst) term in Eq. [ij 
which is diluted by the decay time resolution and mistag 
probability. Systematic uncertainties from these sources 
and from the mixing frequency are absorbed in the sta- 
tistical uncertainties as explained above. Other system- 
atic uncertainties are determined as follows, and added 
in quadrature to give the values shown in Table |l] 

To test our understanding of the decay angle accep- 
tance we compare the rapidity and momentum distribu- 
tions of the kaons and muons of selected candidates 
in data and simulated events. Only in the kaon momen- 
tum distribution do we observe a significant discrepancy. 
We reweight the simulated events to match the data, red- 
erive the acceptance corrections and assign the resulting 
difference in the fit result as a systematic uncertainty. 
This is the dominant contribution to the systematic un- 
certainty on all parameters except Fg. The limited size 
of the simulated event sample leads to a small additional 
uncertainty. The systematic uncertainty due to the back- 
ground decay angle modelling was found to be negligible 
by comparing with a fit where the background was re- 
moved statistically using the sPlot method [16) . 

In the fit each |Ai(0)p is constrained to be greater 
than zero, while their sum is constrained to unity. This 
can result in a bias if one or more of the amplitudes is 
small. This is the case for the S-wave amplitude, which 
is compatible with zero within 3.2 standard deviations. 
The resulting biases on the |Ai(0)p have been determined 
using simulations to be less than 0.010 and are included 




(|)^ [rad] 

FIG. 3. Likelihood confidence regions in the ATa-cps plane. 
The black square and error bar corresponds to the Standard 
Model prediction [Hll]. 

as systematic uncertainties. 

Finally, a systematic uncertainty of 0.008 ps~^ was as- 
signed to the measurement of F^ due to the uncertainty 
in the decay time acceptance parameter /3. Other sys- 
tematic uncertainties, such as those from the momentum 
scale and length scale of the detector, were found to be 
negligible. 

In summary, in a sample of 0.37 fb^^ of pp collisions 
at y/s = 7TeV collected with the LHCb detector we ob- 
serve 8492 ± 97 ^- J/^pK+R- events with K+K' 
invariant mass within ±12 MeV of the (f> mass. With 
these data we perform the most precise measurements 
of (f>s, ATs and Tg in 5° — >■ J/V"/* decays, substantially 
improving upon previous measurements [7] and provid- 
ing the first direct evidence for a non-zero value of AFj. 
Two solutions with equal likelihood are obtained, related 
by the transformation AT.,) i— >■ (tt — 0^, — AF^). The 
solution with positive AFs is 

0s = 0.15 ± 0.18 (stat) ± 0.06 (syst) rad, 
Ts = 0.657 ± 0.009 (stat) ± 0.008 (syst) ps~\ 
AF, = 0.123 ± 0.029 (stat) ± 0.011 (syst) ps~\ 

and is in agreement with the Standard Model predic- 
tion Oil]. Values of (f>s in the range 0.52 < (pg < 2.62 
and —2.93 < (ps < —0.21 are excluded at 95% confi- 
dence level. In a future publication we shall differentiate 
between the two solutions by exploiting the dependence 
of the phase difference between the P-wave and S-wave 
contributions on the K^K~ invariant mass |14) . 
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